Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
内海 隆行*
CIPUS Annual Report, 4 Pages, 1998/00
超高出力超短レーザー光による物質の溶融・蒸発シミュレーションを行うため、多相・多成分の流体が存在する系に適用可能な簡易状態方程式モデルを試作した。これを既存の状態方程式データ(ロス・アラモス研究所のSESAMEとローレンス・リヴァモア研究所のQEOS)と比較することにより、さらに精度を向上するためには膨張領域のコールドEOSの関数型の仮定と、電子EOSのイオン化ポテンシャルの仮定とを詳細化すべきであることを明確にした。
内海 隆行*; 功刀 資彰*; 藤井 貞夫*; 佐々木 明; 赤松 幹夫*
CIPUS Annual Report, 4 Pages, 1998/00
X線レーザー照射によるプラズマ生成過程を連続体運動方程式で記述されるようなマクロスコピックな流体シミュレーションの立場から把握するためには、光と物質その相互作用のモデル化だけでなく、固体ターゲット照射にともなって生じる固・液・気各相の界面での大きな密度変化を精度良く捉えられ、かつ圧縮性高速流れから非圧縮性流れまでの広い範囲の流動現象を統一的に取り扱えることが必要である。矢部によって提案された3次補間擬似粒子法(CIP:Cubic-Interpolated Propagation)は、固・液・気各相を含む圧縮性・非圧縮性流体の統一解法に極めて有効な計算手法と考えられる。筆者らは、CIP法を異なる2種類の物質に対して適用できるように拡張し、多相・多成分の流体が存在する系に適用可能なコードを開発してきた。本報では、そのコードの概要と超短パルスレーザー光を金属に照射する問題に適用した結果について報告する。
内海 隆行*; 功刀 資彰; Zhang, Y.*
CIPUS Annual Report, 0, p.1 - 7, 1996/00
CIP法による双曲型偏微分方程式の数値解析では、移流項計算と非移流項計算が分離して行われるが、移流項計算における数値粘性が小さいことから圧縮性流体運動において現われる衝撃波前後での物理量の不連続を数値的に取り扱うためにVon Neumann-Richtmyer(N-R)型人工粘性項を非移流項計算に導入することが必要となる。これまでCIP法で用いられてきた人工粘性項は、圧縮領域においては衝撃波面における熱エネルギーと運動エネルギーの交換をRankine-Hugoniotの関係を満たすようにし、非圧縮領域においてはゼロとするものであり、空間変数に関しては連続ではあるが微分可能でない関数である。しかし近年CIP法に関連して、微分代数的CIP法、CIP-有限要素、局所補間微分オペレータ法(IDO)のように人工粘性項の微分可能な関数表現を用いることが好ましい数値解析手法が提案されてきている。本報告では、人工粘性項の1次元における解析関数的表現の検討を通して、一般座標・多次元への拡張可能な人工粘性項を提案する。